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ABSTRACT 

Rare earth oxide stabilized zirconias find a wide range of applications. An 
understanding of phase equilibria is essential to all applications. In this 
study, the available phase boundary data and thermodynamic data is 
collected and assessed. Calphad-type databases are developed to 
completely describe the Y 2 03 -Zr 0 2 , Yb 2 03 -Zr 0 2 , and Y 2 03 -Yb 2 03 
systems. The oxide units are treated as components and regular and sub- 
regular solution models are used. The resultant calculated phase diagrams 
show good agreement with the experimental data. Then the binaries are 
combined to form the database for the Y 2 03 -Yb 2 C> 3 -Zr 0 2 psuedo-temary 

INTRODUCTION 

Yttria-stabilized zirconia is one of the most widely used ceramics, with applications 
ranging from thermal barrier coatings to solid oxide fuel cells. There has been interest m 
using other rare earth oxides alone or in combination to stabilize zirconia. These 
different combinations offer improved properties, such as reduced thermal conductivity. 

A survey of phase equilibria data for rare-earth oxide/zirconia systems shows a good 
deal of information on the Y 2 0 3 -Zr0 2 system, but only limited information on the other 
rare-earth oxide/zirconia systems (1-3). Experimental data on two or more rare earth 
oxides in combination with zirconia is even more limited. For this reason, a Calpha 
(CAlculation of PHase Diagrams) assessment is valuable (4). In this method, a series of 
binary databases are combined to form reliable high order databases. 

Due to the widespread use of the Y 2 C>3-Zr0 2 binary system, there are several 
assessments in the literature. Kaufman (5) first approached this system by treating each 
oxide unit as a component. Degtyarev and Voronin (6-7) further simplified this system 



by treating the two cubic phases— cubic Zr0 2 and cubic Y 2 0 3 — as one cubic phase 
with a miscibility gap. Du et al.(l) have thoroughly reviewed all aval able hterature on 
this system through 1990 and developed a Calphad-type database with this miscibi t> 
gap simplification. The group of Yokokawa (8-9) has also calculated the Y 2 0 3 -Zr0 2 
phase diagram and extended these calculations to other rare-earth oxide/zircoma systems. 

In this paper wc rc-asscss the Y 2 0 3 -Zr0 2 system. We do not use the miscibility gap 
simplification, because cubic Zr0 2 and cubic Y 2 0 3 are of different space groups. Also a 
miscibility gap leads to convergence problems from the corresponding multi-component 
systems We also assess the Yb 2 0 3 -Zr0 2 and Y 2 0 3 -Yb 2 0 3 systems and combine these to 
form the Y 2 0 3 -Yb 2 0 3 -Zr0 2 ternary. 


APPROACH 


The first step in the assessment process is to survey all available data m the literature. 
The phases in the Y 2 0 3 -Yb 2 0 3 -Zr0 2 system are given in Table I. 


Table I. Phases in the Y 2 0 3 -Yb 2 0 3 -ZrQ 2 system. 


1 Cl L 1 i X. X. 

Phase 

Abbreviation 

Monoclinic Zr0 2 

mz 

Tetragonal Zr0 2 

tz 

Cubic Zr0 2 

cz 

Liquid Y 2 0 3 -Yb 2 0 3 -Zr0 2 

liq 

Cubic Y 2 0 3 

cy 

Hexagonal Y 2 0 3 

hy 

Zr 3 Y 4 0, 2 

delta y 

Cubic Yb 2 0 3 

cyb 

Hexagonal Yb 2 0 3 

hyb 

Zr 3 Yb 4 0i 2 

delta yb 


As noted, there is a good deal of experimental information on die Y20 3 -Zr0 2 system 
This has been discussed in detail in other publications (1-3) and will only be briefly 
summarized here. The liquidus line of Noguchi et al. (10) was used m this . assessmen. 
There are numerous measurements of the tetragonal/(cubic + tetragonal)/cubic phase 
boundaries (11-17). Several different techniques were used and, m general, these show 
remarkably good agreement. The studies of Stubican et al. (12) are the most complete 
examination^ this system and were used in the assessment. Four invariant points have 
berreported for thfs system-these are listed in Table II. There xs an intermedi^e 
compound of composition Zr 3 Y 4 0 12 , which decomposes to the cubic solid solution of 
ZrOo. A second, pyrochlore-type compound m the Y 2 0 3 -nch portion of the diagram 
been reported by other investigators, but its existence is controversial (18) and it will not 

be considered here. 
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Table II. Invariant points and decomposition temperatures in the Y 2 0 3 -Zr0 2 and Yb 2 0 3 



Temperature 

(K) 

Composition 
x(Y 2 Q 3 ) | 

Reference 

Comments 

Eutectic 
lia = cz + cy 

2688 

0.817 

(12) 


Peritectic 
cv = lia + hy 

2720 

0.778 

(12) 

Introduced, not 
measured 

Euctectoid 
cz — delta v cy 

1650 

0.435 

(13) 


Euctectoid 
tz = mz + cz 

838 

0.039 



(19) 

Difficult to attain 
equilibria at this low 
temperature 

Euctectoid 
cz = delta y + 
mz 

<673 

0.20-0.30 

(12) 

Difficult to attain 
equilibria at this low 
temperature 
Approximate 
euctectoid discussed 
in (12) 

Decomposition 
Temperature of 
delta y 

1655 

0.4 

(20) 


Eutectic 
lia = cz + cvb 

2691 

0.85 

(21) 


Decomposition 
of delta yb 

1898 

0.4 

(20) 



There is very little thermodynamic data for the Y 2 0 3 -Zr0 2 system. Belov and 
Semenov have measured activities in the system at 2773 K (22). For the mtermediate 
compound Zr 3 Y 4 0i 2 , there are some recent first principles calculations of enthalpies at 0 

K (23). 

As noted, there is much less data for the Yb 2 0 3 -Zr0 2 system. The liquidus 
measurement of Rouanet (21) was used in the assessment ^Limited phase ;boundag 
measurements in the Zr0 2 rich side by Stubican et al. (20) were used. The si g 
reported invariant point is given in Table II. This system also contains an intermediate 
compound, Zr 3 Yb 4 0 12 , which decomposes at a higher temperature than the analogous 

Zr 3 Y 4 0] 2 compound (20). 


The third binary needed is the Y 2 0 3 -Yb 2 0 3 system, 
diagram (24) shows a continuous cubic solid solution, 
modeled as an ideal solution. 


The only experimental phase 
Thus this system was simply 
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One study of the Y 2 03-Yb 2 03-Zr02 system was found (25). The major features 
of this are continuous tetragonal Zr0 2 with Y 2 03 and Yb 2 03; continuous cubic Zr0 2 with 
Y 2 0 3 and Yb 2 0 3 ; continuous cubic (Y,Yb) 2 0 3 ; and a continuous line compound 
Zr3(Y,Yb)40i 2 . The presence of the line compound, of course, depends on the 
disordering temperatures in Table II. 


MODELING 

The first step in a Calphad modeling is to determine the lattice stabilities for all 
phases. These were taken from the SGTE pure substance database (26), which m turn is 
from the IVTAN tables (27). These numbers are derived from the phase transition 
temperature and heats. We have modified the melting point of Y 2 0 3 to agree with the 
melting point measured by Stubican (12). The SGTE database includes a sohd/sohd and 
a solid/liquid phase transition for Yb 2 0 3 , but only the temperatures of transition are 
given— no heats are reported. Therefore we used the heats for analogous Y 2 0 3 
transitions. 

Two commercial codes were used for these calculations ThermoCalc and 
WinPhad**. Both express the solutions as: 

G m (x,T) = X x i° G > (T) + RT£x ; InXj + ~G m 0 ) 

i=1 

ex G m = x 1 x 2 ^ j L(x,-x 2 ) J ( 2 ) 


The first term in equation (1) is the sum of the Gibbs energies of each component, the 
second term is the ideal entropy of mixing, and the third term is the excess Gibbs energy 
Both codes express the excess Gibbs energy as a Redlich-Kister polynomial (28). In all 
cases, only the first one or two terms were needed, which corresponds to the regular and 
sub-regular solution models, respectively. 

The solution parameters, were derived from the selected experimental data using the 
optimizing modules of ThermoCalc and WinPhad. The most prominent features of the 
phase diagrams were first developed and then the fine features added. The goal is to 
develop a description of the entire system with a minimum number of parameters. As 
noted, binaries can be put together to form ternaries. In general a limited number of 
ternary interaction coefficients are necessary. 

RESULTS 

Figure 1 is the calculated Y 2 0 3 -Zr0 2 psuedo-binary. In general agreement with 
experimental data is good. The four invariant points are reproduced as is the disordering 
temperature for Zr 3 Y 4 0i 2 . The (tetragonal + cubic)/cubic phase boundary is x(Y 2 0 3 ) 
-0.02 high at the lower temperatures. A better model should help reproduce this feature. 
Current work focuses on the modeling of this system with the ionic model (29). 


' ThermoCalc Software, Stockholm, Sweden 
*’ CompuTherm LLC, Madison, Wisconsin 
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Figure 2 is the calculated Yb 2 0 3 -Zr0 2 psuedo-binary. Agreement with the limited 
experimental data is good. Note the intermediate compound decomposes at a hig er 
temperature than that in the Y 2 0 3 -Zr0 2 system. 



* Noguchi 

0 Stubican (1978) 
x Ruhie 

* Srivastava 

A Stubican (1984) 

* Yoshikawa 

* Pascual 


Figure 1. Calculated Y 2 0 3 -Zr0 2 phase diagram. The data points are the experimental 
points. 



# Rouanet 
D Corman 


Figure 2. Calculated Yb 2 0 3 -Zr0 2 psuedo-binary phase diagram. Data points are 
experimental data. 
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Figure 3. Calculated Y 2 03 -Yb 2 03 -Zr 02 isothermal section at 1673 K. 

Figure 3 is the calculated psuedo-temary section at 1673 K. The major experimental 
features of this ternary (24) have been reproduced. There are continuous tetragonal Z 2 , 
cubic Zr0 2 , and cubic (Y,Yb) 2 0 3 phases. Note also the intermediate line compound 
Zr 3 (Yh,Y) 40 i 2 , does not extend to the Zr 3 Y 3 0i 2 side as this compoun as area y 
decomposed at 1673 K (see Figure 1). 

SUMMARY AND CONCLUSIONS 


The Y-i0 3 -ZrOi, Yb 2 0 3 -Zr0 2 , and Y 2 0 3 -Yb 2 0 3 psuedo binaries have been assessed. 
Av^abk phase boundaiy and thermodynamic has been used to derive the soludon 
parameters P for a Calphad-type database. The oxide units are used as components, 
resultant phase diagrams show good agreement with the experiment daUu Th«m *e 
three binaries have been combined to form a database for the Y 2 3 23 P 

ternary. 
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